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Abstract. Transmission of light through an atomic sample placed between 
crossed polarizers in a magnetic field is jointly determined by Faraday rotation 
and net absorption: transmission increases with rotation and decreases with 
absorption. Both rotation and absorption are proportional to the atomic den- 
sity No, hence, in a certain range of No, the two effects may compete yielding 
a distinct density dependence of the transmitted light. We have studied such 
dependence in rubidium vapor for No ~ 6.0 X 10 9 — 3.0 X 10 11 at/cm 3 with 
resonant laser light. We present interpretation of the competition effect and 
discuss its possible application for atomic density determination. 



1. Introduction 

Studies of magneto-optical effects in atoms and molecules have a long tradition. 
The advent of tunable lasers stimulated renewed interest in the nonlinear magneto- 
optics. Numerous papers are describing various aspects of this nonlinearity (see 
referencepQ for recent review). In the systems with ground state (or other long- 
lived lower levels) of angular momentum J ^ 0, the magneto-optical nonlincaritics 
are due to optical pumping and Zeeman coherences. Optical pumping changes the 
saturation conditions for a given transition by affecting the sublevel populations 
whereas the Zeeman coherence allows interference of quantum paths linking differ- 
ent sublevels. These effects give rise to an optical anisotropy which can be easily 
studied experimentally 

One important application of magneto-optical studies is the determination of 
optical density and thus a measurement of oscillator strength (dipole moment) of a 
given transition when the number density is known or, inversely, of atomic densities 
when the transition probability is known. Examples of such applications of linear 
Faraday effect are the measurements of relative oscillator strengths in neon |2j and 
sodium [3] and of atomic vapor density of optically thick potassium vapors @|. 

The present work describes experimental and theoretical analysis of transmission 
of resonance laser light by atomic vapors positioned between crossed polarizers in 
a weak longitudinal magnetic field (forward scattering - FS) under conditions of 
nonlinear light-matter interaction. In general, the FS signal, i.e. the transmitted 
light intensity vs. magnetic field, depends on the polarization anisotropy both 
in atomic absorption and dispersion. In such context FS has been analyzed for 
conventional light sources in early works of G. W. Series and coworkers and for 
laser light by W. Gawlik et al. fi. . For a medium of length L, placed between crossed 
polarizers and characterized by refractive indices n± and absorption coefficients k± 
for the (T^~ polarized light, the intensity of the forward-scattered monochromatic 
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light of wavenumber k is given byPJ 

(1.1) I FS = ^[e- K + kL - e-"- kL f + Jo sin 2 [ " + ~ — kL} e -^++^ kL , 

where Iq is the incident light intensity. The FS signal is non-zero only if the medium 
possesses some optical anisotropy between the a + and <j~ light polarizations. An 
important cause of such anisotropy is a longitudinal magnetic field which splits 
the n±(u>) and k±(lu) profiles. The first part of Eq. (|1.1|) describes the effect of 
medium dichroism whereas the second part is due to its birefringence weighted by 
a net absorption. For the light frequency to tuned exactly to the atomic transition 
frequency u>q and for a linear splitting of atomic energy levels due to magnetic field, 
the appearance of a nonzero FS signal is caused exclusively by the second part 
of expression (|l.ll) . i.e. by the birefringence. The exact shape of the FS signal 
dependence on the magnetic field is affected by a combined effect of birefringence 
and total absorption. In the case of nonlinear interaction with the light beam, 
this dependence may have the form of a narrow sub-Doppler, or even subnatu- 
ral, resonances (see, e.g. references The two properties of the medium, 

birefringence and total absorption, scale differently with the atomic density which 
results in specific density/temperature characteristics of the FS signals discussed 
below. We hope that these characteristic density dependences would be useful for 
determination of the optical density of atomic samples. 

The detailed mechanism of nonlinear FS signals has been theoretically analyzed 
in several papers (see, e.g. references[7| |S1 EB ^2 El) at different levels of 
accuracy. Below, we present a nonperturbative treatment that for low J values 
yields analytical formulae. 

In the next section we describe our setup and present experimental results. In 
Sec. 3 we present theoretical calculations and interpretation of the experiment. 
Finally, we conclude in Sec. 4. 

2. Experiment 

2.1. Apparatus. The experimental setup is shown in Fig. 1. We used an external 
cavity diode laser tuned to the hfs component F g = 3 — > F e — 2 of the 85 Rb D-y 
line. The laser frequency was stabilized by an external stabilization system based 
on Doppler-free dichroism |13|. The light beam of 2 mm diameter was linearly 
polarized and its power was varied by rotation of a A/2-waveplate in the range to 
4 mW. The laser light was passing through a 5 cm long cylindrical glass cell of 2 cm 
diameter filled with the natural mixture of 85 Rb and 87 Rb and no buffer gas. The 
cell was heated by a nonmagnetic heater to temperatures between 20 °C and 60 °C 
which yielded saturated vapors with densities in the range 6.0 x 10 9 -j- 3.0 x 10 11 
at/cm 3 . A pair of Hclmholtz coils produced the longitudinal magnetic field of -3 
G to 3 G value. The cell and the coil system were surrounded by a three-layer 
magnetic shield made of mu-metal. As assessed from the resonance line-width, the 
shielding efficiency was better than 10 4 . Polarizer P and crossed analyzer A were 
mounted outside the shield|S]- In that geometry, detector D placed behind the 
analyzer registered only the FS signal due to rotation of the polarization plane of 
the transmitted light and/or its ellipticity. Part of the incident beam was directed 
to the saturation spectroscopy system which served as frequency reference. 
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2.2. Results. The FS signals, Ifs(B), were recorded as a function of the magnetic 
field B for various temperatures T and light intensities Iq. The narrowest features in 
signals are caused by nonlinear Faraday effect due to ground-state optical pumping 
and Zeeman coherences. A typical plot of that part of the FS signal is presented 
in Fig. 2a. The amplitude of the signal, defined as Aps in Fig. 2a, depends on 
the light intensity and concentration of the atomic vapors. In Fig. 3 we summarize 
the results of the amplitude vs. temperature T dependences taken for various light 
intensities Iq. 

For the studied range of intensities between 1 mW/cm 2 and 6 mW/cm 2 , the 
signals exhibit characteristic temperature dependence: for low temperature their 
amplitudes increase with T then reach their maxima and subsequently decrease with 
T . The temperature of the maximum amplitude shifts with the light intensity. For 
high temperatures almost no light is transmitted through the cell between crossed 
polarizers. For very high light intensity (above 100 mW/cm 2 ), the maximum shifts 
even beyond the range of temperatures available with our setup (Fig. 4). 

We have verified that the amplitude of the FS signal is a nonlinear function of the 
light intensity. This could be expected as the intensities we were using were close 
to, or higher than, the saturation intensity which for rubidium is 1.6 mW/cm 2 . It 
can be thus expected that the interplay between rotation and absorption is affected 
by different saturation behavior of the dispersive and absorptive atomic properties. 
Such different saturation effects are manifest in Fig. 3 by a shift of the maxima 
towards higher temperatures with increasing Iq. 

3. Interpretation 

The properties of the maxima in the Aps(T) dependences, described above, 
are independent of the detailed mechanisms responsible for the FS signals. The 
maximum occurs for linear, as well as for nonlinear Faraday rotations. For nonlinear 
Faraday effect, it appears no matter if the FS signals have natural or sub-natural 
widths, i.e. if they are due to population redistribution or Zeeman coherences. The 
detailed nature of atomic birefringence and absorption mechanism affects only the 
range of the concentrations (temperatures) where the maximum is visible. 

This property can be illustrated by the following simple arguments. For exact 
resonance, S = 0, only the second part of Ea. (|l.l(l differs from zero, hence 

(3.1) ^ = sin 2 (aiVo) exp(-/3A/ ), 

Jo 

where a/(3 — l/2(n + — ti_)/(k + + k_) and No is the atomic concentration. For 
a given atomic transition, cell length, magnetic field and light intensities, a and 
(3 are constants mutually related by the Kramers-Kronig relation [T3|. Ips/Io is 
thus an universal function of Nq, parameterized by the ratio a//3 but independent 
of the detailed form of a and (3. Generally, the specific value of the ratio depends 
on particular dependences of a and f3 on atomic energy level structure, mechanism 
of interaction with the external fields and their intensities. However, for the value 
of the magnetic field corresponding to maximum rotation Aps the splitting of 
atomic magnetic sublevels is nearly equal to their width which results in a/ [3 of 
the order of 1. By virtue of the Kramers-Kronig relation, the value of a/ (3 at 
a given B is nearly insensitive to the details of the atomic structure and to the 
nature of optical anisotropy. With a/ (3 = 1, expression (|3.1|) represents a strongly 
damped oscillatory dependence on No, dominated by the first oscillation. These 
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oscillations arc nonlinear counterparts of those described in reference [2j- Note, 
that the appearance of the characteristic maximum of the Afs dependence on Nq 
is possible only with the FS methodology and cannot be observed by measuring 
only the rotation. 

For more quantitative interpretation of our experimental results we use semiclas- 
sical approach of forward scattering and density matrix formalism ^2] to calculate 
the refractive indices n and absorption coefficients k 

(2a) ,-l=^y^H 



(2b) K= ^^M 

1 ' Ee 

where {pflh)v are the velocity averaged optical coherences for the i-th transition of 
a given polarization, 5ft and 9 denote their real and imaginary parts, respectively, 
di is the dipole moment of the particular transition, E the amplitude of the light 
electric field, and eo is the electric permittivity of free space. 

In this paper we consider the simple model of atomic structure J g = 1/2 — > J e = 
1/2 with two sublevels in the ground state and two sublevels in the excited state 
the, so called, X model. The er + component of the propagating beam excites atoms 
from the ground-state sublevel \a) (m = —1/2) to the excited-state sublevel \B) 
(to = 1/2). Analogically, the a~ component causes transition from |6) (to = 1/2) 
to |^4) (to = —1/2). The great advantage of this simple model is its analytical 
solvability. Despite its simplicity, the X model gives good qualitative agreement 
with the experimental results. Furthermore it reproduces well the observed depen- 
dence of the FS amplitude on the properties of the atomic medium, in particular on 
the dipole moment and concentration of atoms in the ground sublevels in thermal 
equilibrium. 

To check to what extent the density dependence of Aps(No) is insensitive to 
particular model, and to verify that the X model does not oversimplify the in- 
terpretation, we have additionally applied another model using the A-type atomic 
structure (Fig. 5b). The second model allows proper analysis of the coherence 
effects, such as coherent population trapping (CPT) and their consequences 
for magneto-optical properties of the medium. However, analytical solutions and 
velocity integration are much more difficult within the second model. For this rea- 
son, in this paper we concentrate exclusively on the X model. Nevertheless, we have 
verified by numerical calculation with the A model that Zeeman coherence effects 
do not change qualitative predictions of the X model concerning the temperature 
dependence of ApS(N). This is due to the fact that the IfS(B) signals simulated 
by both models are qualitatively very similar to the experimental ones, as seen in 
Fig. 2, so the a and (5 quantities in Eq. (2) remain very close, no matter what 
detailed model is being used. 

Within our model, the refractive indices and absorption coefficients for both 
circular polarizations are calculated using equations 1 2 ;i I and (|2b|) as 

/o \ i dft((p aB ) v ) d5t((p bA )v) 
(3a) n + - 1 = , n- - 1 = , 
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(3b) K+ = Ee ' K = Ee ' 

The semiclassical Hamiltonian of the system is given by 

(4) ^ 

H = huJo(p AA +pBB) + -^9g^B(Paa-Pbb) + -^9e^B(PAA-pBB) + dE(p aB +PbA + h.C.), 

where pa denotes population of the i-th state, g g , g e are the Lande factors of the 
ground and the excited states, respectively, ujq is the transition frequency in zero 
magnetic field, and u>b is the Larmor frequency lob = PbB/H with ps being Bohr 
magneton. 

Evolution of the density matrix is described by the master equation 

(5) p = - l -[H,p]-tp, 

where Tp stands for the relaxation operator. 

Formula JHJ) yields the following equations describing evolution of the optical 
coherences in the X model 

(6a) p aB =i{U()+Uj'B+ ft) PaB + i(3(paa ~ PBb), 



(6b) p bA = i(u -lo' b + i^)pbA + if3(pbb - Paa), 

where ui' B = \ (g g + g e )u B is the effective magnetic shift of the resonance frequency, 
/3 = Ed/h is the Rabi frequency and T is the relaxation constant of the optical 
coherences. 

Further calculations are performed in the rotating-wave and steady state approx- 
imations. Under such conditions the optical coherences become 

(7a) PaB = -T- , , . r (Paa - PBB), 

(7b) PbA = -- x f . .A Pbb ~ Paa)-, 

o — uj b + 11 

where p a B and p^A are slowly-varying coherences in the rotating-wave approxi- 
mation and 8 = luq — u>. Although only p a B is calculated below, the analogical 
calculation can be performed for ps a - 

Relation (|7a|l couples the optical coherence with atomic populations. Using Eq. 
JSJl to calculate atomic populations and inserting them to (|7a[l one can write 

te\ - a fx > -"p\ M A (S) + 2rG 

(8) PaB - f3n eg (6 + u B - iT) Ma{s)Mb{S) + rG[MA{5) + Mb{5)) + 2rG 2 ■ 

where n eq is the thermal equilibrium atomic population of the ground state and 
other symbols are defined as 

M A {S) = (6-uj'b) 2 +T 2 , 

m b (s) = (s + aj' B ) 2 + r 2 , 

7e 
3 7g 
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with 7 g and 7e being the relaxation constants of the ground and excited states, 
respectively. Relation (jHJ is an analitical expression valid for any light intensities 
and allowing easy analysis of the dispersive and absorptive responses of our model 
system to light and magnetic field. Perturbative solutions for optical coherences 
can be obtained in the low light power regime, G <C 1, by expanding relation (|SJ) 
into power series of G. 

Relation (jHJl contains the ratio of two polynomials of S (they also depend on 
lo' b ). At this point atomic movement needs to be taken into consideration, i.e. S 
has to be corrected for Dopplcr shifts (5 — > S — kv) and expression JSJ needs to be 
integrated over the Maxwell distribution of the velocities v. Generally, this is the 
place when one sacrifices analyticity of solutions for a single velocity class. However, 
the benefits of analytic expressions can be retained by expanding © into partial 
fractions 

[1 PaB ~ neqn S + 5 1 (co> B ) + 6-5 1 (lo' b ) + S + 6 2 (lo< b ) + 5-6 2 (uj> b ) U 



where 5 2 12 = B 2 - 7 2 - rG ± ^Jr 2 G 2 - 4B 2 ( 7 2 + rG) - 2rG 2 are the roots of the 
denominator of Eq. JSJ, while A±, B± are the expansion coefficients equal to 

1 6f T aS'l + bdi =F c 1 (5| =p + W 2 T c 

±= 2 ^(<5 2 -<5 2 2 ) ' ± = _ 2 " <5 2 (5 2 -<5 2 ) ' 

where a = B + 17, 6 = 2rG - (a*) 2 , and c = -a*(2rG + |a| 2 ). One thus obtains 

n eq (3N , f°° e?q>(-v 2 /u 2 ) dv 

{PaB)v = 7= ^+(W B J 

(10) 



-00 + 
exp(-w 2 /u 2 ) 
x 5 - kv - 5\{ui' B ) 
a exp(— v 2 /u 2 ) dv 
x S - kv + 5 2 (uj' b ) 
exp(— v 2 /u 2 ) dv 
5 - kv - 8 2 {uj b ) 

where u is the average atomic speed. Writing down the coherence in form <|1U|) 
allows a clear interpretation of each contribution. The functions A±(oj' b ), B±(u>' B ) 
are responsible for narrow structures in the FS signals and describe the nonlinear 
Faraday effect, while the integrals describe wide pedestals of Ipg{u:' B ), barely visible 
in a narrow scan of the magnetic field. For calculations of the integrals in (|10|l . 
standard tabularized functions, the plasma dispersion function |17| or the complex 
error function |18| . can be used. 

As it was mentioned above, the FS signal registered with the crossed polarizers 
arrangement results from competition between rotation of the polarization plane 
and net absorption. The competition is seen in Ea. (|l.l|) which for 5 — 0, takes the 
form 

(11) I FS = i sin 2 ( n+ ~ 2 n ~ kL) exp(~2nkL), 

where k = k+(cjb,(5 = 0) = k-(ojb,S = 0). Relation (|10|) was used for simulation 
of the FS signals with parameters corresponding to the experimental conditions: 
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d = 2.25 xl0~ 29 Cm, u = 200 m/s, g e = g g = 1, ~/ g = 2tt x 0.1 MHz, 7e = 
27r x 5.74 MHz, and T = 2tt x 2.92 MHz. The ground-state relaxation rate 7 S 
reflects finite atomic transit time across the laser beam. In the D\ line of rubidium, 
there are two transitions within the Doppler width that can be excited by light: 
F g = 3 — > F e = 2 and F g = 3 — > -F e = 3 but since in our experiment the laser 
was tuned only to the F g = 3 — > F e = 2 component, the contribution of another 
transition has been neglected. Apart from a very simplified atomic structure, the 
present analysis also does not take into account specific relaxation mechanisms, e.g. 
velocity-changing collisions^] and/or radiation trapping |20I which may affect the 
specific ratio a/ (3 in Eq. (|3.1() . 

Results of the simulations with our X model of the FS signal amplitude vs. 
atomic concentration are shown in Fig. 6 for various light intensities. 

As the experimental results were recorded vs. vapor temperature, the atomic 
density in Eq. I|10|) has to be expressed as a function of the vapor temperature which 
is possible if the temperature dependence of rubidium vapor pressure is known. 
The issue of vapor pressure is crucial for many experimental situations, yet exact 
measurements of vapor pressures are difficult and subject to various systematic 
errors. A common procedure is to use one of the existing compilations which fit 
experimental data obtained in various temperatures by simple analytical formulae 
that can be used for extra- or interpolation to other temperatures. In principal, 
such fits are based on the statistical physics and chemistry laws (Arrhenius law, 
ideal gas law) with possible corrections for departures from the ideal gas model 
and/or for specific metal evaporation conditions. Particular care should be taken 
close to phase-transition temperatures, as e.g. near the melting point. This is the 
case of our experimental conditions as the melting point of rubidium 39.3 °C occurs 
just in the middle of our temperature range. 

Assuming ideal gas and starting from the Arrhenius relation, the saturated vapor 
pressure could be determined by formula 

(12) )ag w p(T) = A-~. 

Such a dependence is used, e.g. in reference [2~T]. which for rubidium gives A = 
4.875 and B = -4215 below the melting point and A = 4.312 and B = -4040 
above the melting temperature. Another widely cited compilation |22j applies extra 
corrections terms and uses relation 

(13) log 10 p(T) = A' + ^ + C'T + D' log 10 T, 

where for rubidium A = 94.05, B' = 1961, C = -0.03772, and D' = 42.57 below 
the melting temperature and A 1 = 15.88, B' = 4330, C = 0.00059, and D' = -2.99 
above the melting point. From Eqs. l(T2*j l and iJTSjl . with the help of ideal gas equa- 
tion, pV = N a tkBT (N at is the number of atoms the gas, V is its volume, and ks 
the Boltzmann constant), one can obtain the final concentration-temperature rela- 
tion Nq(T). Reliability of such density determination depends obviously on quality 
of experimental data used for a given analytic fit. Reference j23| presents a thorough 
comparison of several measurements of rubidium vapor pressure vs. temperature 
which appears to be fairly consistent with the Nesmeyanov's prediction[22] • On the 
other hand, the CRC compilation 21 is based on not too many experiments and, 
to our best knowledge, has not been verified in a wide range of temperatures. 
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Our experiment interpreted in terms of the X model yields for T « 50 °C the 
density that is consistent with the data given by|5J and by 25% higher from the 
value given in reference |22j. We hope that after more detailed interpretation, the 
described method could allow accurate measurements of the atomic vapor densities 
and possibly improve existing literature inconsistencies. 

Comparison of Fig. 3 with Fig. 7 shows that while there is not a perfect agree- 
ment of the simulated signals with the experimental ones, they are very similar and 
exhibit the same general character. For low temperatures the FS amplitude scales 
as a square function of the concentration No. This part of the dependences is de- 
termined only by the dispersive properties of the medium. For higher temperatures 
the absorption starts to play an important role and cannot be neglected anymore. 
The interplay between dispersion and absorption produces the maximum of the FS 
amplitude signal. For even higher temperatures, the absorption of the transmitted 
light becomes dominant and the amplitudes of the FS signals decrease exponentially 
with concentration Nq. In addition to this temperature/density dependence, the 
FS amplitude rises and shifts with increase of the incident light intensity. The shifts 
of the maxima positions for various light intensities are due to different saturation 
behavior of the dispersion and absorption coefficients which affects the a/ (3 ratio 
in Eq. (|3.1I) . This is a very interesting effect that deserves independent systematic 
study, preferably with more complex models of atomic structure. 

4. Conclusions 

The specific effect of competition between polarization rotation and absorption, 
described above, shows up in the FS signals. The FS intensity, given by Eq. (JTDJ, 
increases with the medium birefringence (as sin 2 of the Faraday angle) and decreases 
exponentially with the medium absorption. Competition of these opposite trends 
is responsible for the observed maxima of the dependences seen in Figs. 3 and 
7. Appearance of this competition is an universal feature, independent of whether 
the light-matter interaction has a linear or a nonlinear character. The possible 
nonlinearity only affects the absolute scale of the rotation and absorption and, 
consequently, the range of densities where the competition occurs. For alkali metal 
vapors and for light intensities of the order of 1 mW/cm 2 , the resonant nonlinear 
Faraday effect allows observation of such a competition already in a low-density 
range, about Nq ~ 10 10 at/cm 3 , which requires temperatures not much higher 
than the room temperature. The onset of optical nonlinearity not only lowers the 
density at which the characteristic maximum of the Afs(Nq) dependence appears 
but results in its shift with the light intensity. This is due to different saturation 
behavior of the dispersive and absorptive atomic properties. 

We note that the main features of the discussed dependences do not depend on 
the detailed mechanism of the nonlinearity. We have considered two physically dif- 
ferent models in which the nonlinearities are caused by two different mechanisms: 
either they are due to redistribution of populations between Zeeman sublevels (op- 
tical pumping), or are due to the coherence established between the sublevels. Both 
models gave very similar results and reproduced our experimental observations. The 
very distinct result of the competition between absorption and rotation, i.e. the 
maximum of the forward scattering amplitude signal, appears at well determined 
value of the atomic density. We believe that existence of such characteristic point 
whose position is well defined function of the atomic properties of the medium and 
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scattered light intensity should allow precise determination of the atomic number 
densities or transition probabilities. 
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Figure 1. Experimental setup. BS is the beam-splitter, P, A are 
polarizer and analyzer, respectively, A/2 is A/2-waveplate and D is 
a photodetector. 



12 



SZYMON PUSTELNY, WOJCIECH LEWOCZKO, WOJCIECH GAWLIK 




3 

CD 



J 0.4 




Figure 2. Typical forward scattering signal, Ifs(B), recorded for 
cell temperature 29.5 °C and light intensity Iq = 2.2 mW/cm 2 (a). 
FS signals calculated with the X model (solid line) and A model 
(dashed line) (b). Aps denotes the amplitude of the FS signal. The 
theoretical signals correspond to incident light intensity Iq = 0.6 
mW/cm 2 . They are presented in domain of Larmor frequency over 
effective width (r e / = T for the X model and T e f = j g for the A 
model) and are normalized to the same Aps- 
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Figure 3. The amplitudes of the FS signals (Aps) vs. tempera- 
ture for various intensities of the incident light. 
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Figure 4. The amplitude of the FS signal (Aps) vs. temperature 
for very high light intensity, Iq — 130 mW/cm 2 . 
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Figure 5. Energy level structure for a) the X-model, b) the A- model. 
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Figure 6. Results of numerical simulations of the FS signal am- 
plitudes vs. atomic concentration for the same incident light in- 
tensities as used in experiment. 
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Figure 7. Results of numerical simulations of the FS signal am- 
plitudes vs. vapor temperature for the light intensities as used in 
the experiment (compare with Fig. 3). 
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